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We study the effect of extrinsic noise in metabolic networks. We introduce external random 
fluctuations at the kinetic level, and show how these lead to a stochastic generalization of standard 
Metabolic Control Analysis. While Summation and Connectivity Theorems hold true in presence 
of extrinsic noise, control coefflcients incorporate its effect through an explicit dependency on the 
noise intensity. New elasticities and response coefficients are also defined. Accordingly, the concept 
. of control by noise is introduced as a way of tuning the systemic behaviour of metabolisms. We 

' argue that this framework holds for intrinsic noise too, when time-scale separation is present in the 

system. 

1. INTRODUCTION 

Stochastic fluctuations represent an important contribution to complex behaviours of biological systems. Stochas- 
' ticity appears as a fundamental dynamical mechanism, which does not only generate phenotypic diversity Q, but also 
^0 . plays a major role at many different levels. Selection of alternative pathways in epigenetic switches 0], or synchroni- 
Q ' sation of multicellular systems (33| . are just two of the many processes which not only are influenced by stochasticity, 
■ . but seem to use it to perform in optimal way [sS] . 
1^ ' In molecular biology two classes of stochastic fluctuations are particularly relevant. The flrst class is related to 
the low copy number of chemical species. In particular, if N is the number of molecules in the system, fluctuations 
in N lead to an associated statistical noise with intensity of the order of N~^/'^. While continuous deterministic 
y—{ . descriptions in terms of average concentrations can effectively capture_the relevant dynamics for N large, when N is 
^ ' small fluctuations may become huge, and noise cannot be neglected .9']. Gene regulation is a typical example, as it 
. may be affected by large fluctuations due to the low copy number of transcription factors . Fluctuations associated 
with the intrinsic discreteness of the collisional processes among single molecules are usually referred to as intrinsic 
(or internal) noise. 

On the other hand, the behaviour of a biochemical system depends also on a number of control parameters. Some 
of these relate directly to the macroscopic environment, such as for instance illumination conditions, or pH levels. 
Others have a more microscopic origin, but still exhert a control on the sytem which is independent of its intrinsic 
: dynamics. In gene regulatory networks, for instance, factors acting globally on all genes, such as abundance of 
^ ' RNA polymerase, can change the global efficiency of transcription factors, and ultimately contribute to tuning gene 
expression levels [26| . Similarly, reaction constants, enzyme activities, or input signals critically control the functional 
behaviour of metabolic pathways. In thise sense, all these parameters are external to the system. Fluctuations of 
external parameters define a second type of noise, usually referred to as extrinsic, or external. 

The effect of extrinsic noise can be highly non-trivial and counterintuitive. Extrinsic noise in homogeneous chemical 
systems has been shown to provoke noise-induced transitions [22l| , and recently it has been identified as a mechanism 
for creating and sustaining spatio-temporal patterns in spatially extended systems . In gene networks the interplay 
between intrinsic and extrinsic noise has been recently analysed in [2^, [2^ [2^ . 

In this paper we plan to extend these findings to metabolic networks. We shall focus on the case when the copy 
number of molecules is large enough for intrinsic noise to be negligible, so that a continuous description in terms 
of average concentrations is feasible, but at the same time the system experiences extrinsic stochastic fiuctuations. 
Deterministic rate equations become thus stochastic differential equations. 

Given the description in terms of rate equations, it is ideally desirable to define other, more "systemic" approaches, 
which explain how global properties of the pathway, such as fluxes and concentrations, depend on local variables, such 
as enzyme activities. To this aim, in the specific case of metabolic networks a useful strategy -the so-called Metabolic 
Control Analysis (MCA) [iBj- has been developed, and has become nowadays a popular quantitative framework for 
investigating control and regulation of metabolisms. 

Mathematically MCA is a sensitivity analysis, and is based on a proper manipulation of the rate equations. It 
consists in perturbing the parameters of the system, and evaluating the corresponding change in steady state fluxes 
and concentrations. Perturbations need to be small so that a linear approximation can hold, even if in principle the 
analysis can be carried on including higher order terms. The effect of the perturbation is represented by a control 
coefficient, which expresses the control that a certain rate exerts on the steady state variable of interest. 
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A fundamental property of control coefficients is their interdependence. This is formalized in the so-called Sum- 
mation Theorems 15] for both flux and concentration control coefficients: For any flux, the sum of the flux control 
coefficients of all the rates in the pathway equals one. Similarly, for any metabolite, the sum of the concentration 
control coefficients of all the rates in the pathway equals zero. These properties demonstrate the actual complexity of 
the system. In fact, control on systemic variables is typically distributed among all the enzymes in the system, being 
the respective control coefficients often of the same order of magnitude and constrained by the Summation Theorems. 
This implies that none of them alone is capable of determining the systemic properties of the network [sj • 

Control coefficients are also responsible for mediating the way local perturbations climb up to the systemic level. 
This is the content of the so-called Connectivity Theorems and Partitioned Response relations. The response of the 
system is the sum of the product of the local sensitivity of the individual rates to perturbations (elasticities) times the 
control that those rates exert on the system variable of interest. Connectivity Theorems and Partitioned Response, 
together with the Summation Theorems, define the mathematical structure of Metabolic Control Analysis. 

We present here a novel framework, which allows the description of extrinsic fluctuations within Metabolic Control 
Analysis. Our approach shows that when considering extrinsic stochasticity, Summation Theorems are still valid 
in the standard form, but with control coefficients explicitly dependent on the noise intensity. The same is true 
for Connectivity Theorems and Partitioned Response, and for elasticities and response coefficients. Our findings 
suggest a reinterpretation of noise as an "active" player in the global dynamics of the network, rather than as a mere 
uncontrollable variable. On the one hand, noise may allow the metabolic pathway to access dynamical regions which 
would be otherwise out of reach if the dynamics were purely deterministic. On the other hand, it may provide us 
with an alternative mechanism to exert external control on the regulation of metabolisms. 

In Section 2 we define the stochastic rate equations associated with extrinsic noise, and we illustrate how the 
corresponding stationary probability distribution shows non-trivial modifications with respect to the deterministic 
case. Section 3 is the core of the paper. When extrinsic stochasticity is implemented, the general principles of MCA 
still hold, and incorporate in a natural way the effect of noise. An explicit illustration of our approach is presented in 
Section 4, and we finally conclude in Section 5. 



2. KINETICS WITH EXTRINSIC STOCHASTICITY 



Consider a metabolic network with N metabolites and M reactions. In the homogeneous (well-stirred) situation, 
spatial dependencies are negligible, and the network is dynamically represented in terms of the following set of ordinary 
differential equations: 



dc 
'dt 



f(c,/x), (1) 



with 



M 



/, = z = l,...,7V. (2) 

The N X M matrix Sij is the so-called stoichiometry matrix, and Vj = Vj (c, fi) are reaction rates. Here c = {ci, cn} 
represents the vector of concentrations, and fi — {^i, ^_r} is a vector of R parameters, such as reaction constants, 
enzyme concentrations, illumination conditions, input signals, etc. 

We are interested in the situation when these parameters behave dynamically, and in particular exhibit random 
fluctuations. Experiments can be designed on purpose to investigate the effect of extrinsic fluctuations. For instance, 
in [l^ the effect of stochastic fluctuations on the substrate injection rate in a model for glycolytic oscillations is 
analyzed, and a new noise-dependent dynamical regime of oscillations is identified. On the other hand, extrinsic noise 
can originate from intrinsic noise after some reduction procedure has been carried out. For instance the derivation of 



the Michaelis-Menten kinetics relies on the assumption of rapid equilibration of the enzyme |15l.l20|. As a consequence, 
the dependency on the enzyme concentration is parametric in the Michaelis-Menten form, and intrinsic fluctuations 
on enzyme activity are then perceived as extrinsic in the reduced model. Enzyme activity can fluctuate because of 
the stochastic nature of gene expression ^23| , or due to protein conformational changes 10 1 . 



Hence, we model extrinsic random fluctuations on the parameters fj, by setting for i = 1, i?: 

f^i ^ fJ.i{t) ^ fi.i+^i{t), (3) 
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where Hi is the average value of Hi, and the stochastic process (,i{t) satisfies the equation 

5 = + (4) 
at Ti Ti 

Eq. ^ defines the well-known colored Ornstein-Uhlenbeck noise [l^]- The parameters play the role of the inverse 
correlation time of the noise, while the parameters £i are the noise intensities. These account for different parameters 
to respond with different sensitivities to the external fluctuations. The case when one or more parameters do not 
fluctuate at all, is selected by choosing the appropriate e^'s equal to zero. Furthermore, the assumption that 

£4(^1)^0 as fii for i — l,...,R (5) 

guarantees that fii stays positive, and centered around /Ij . Finally the noises Ci (t) a-re assumed to be Gaussian, with 
zero average, and correlations given by 

m)CAi')) = ^s^A^ - ^'), hj = 1, -^R- (6) 

By performing a linearization of the function f (c, fi) around p. — {/ii, P.r}, Eq. ([T|) is replaced by 

R 

= fJr fA -4- ' 

dt 



= /,(c,/i) +^5,j(c,/i)Cj(t), i = l,...,iV (7) 



3 = i 

with 

9i] c, fi) = — 



(8) 



Eq. ([9]) represents a stochastic differential equation driven by multiplicative Ornstein-Uhlenbeck noise. For the sake 
of simplicity we make the further approximation that the noise correlation times Ti 's be much smaller than any other 
time scale in the system. In this case the noises can be considered as virtually white. The process originally deflned 
by ID), ([3]), and ^ becomes then to the first order in Si and zero order in t^: 



dt 

with 



R 

5 = /,(€,/!)+ ^£yV,(c,/i)0 (i), i = l,...,N (9) 



{Um{t'))^2S,,5{t~t'), i,j^l,...,R. (10) 



Assuming the noise as virtually white is legitimate in all those cases when fluctuations are rapid, so that the zero Ti 
expansion is justified. This may be the case in many different experimental situations, as it is highlighted for instance 
in ^IS*!. Also enz yme activities fluctuations as related to conformational changes have been proven recently to have 
fast components |lOl[l7]. Environmental fluctuations, such as light fluctuations, can be experimentally realized so as 
to be very rapid as well. On the contrary, if fluctuations are not faster than the other typical processes of the system, 
the present scheme breaks down. In this case higher order terms in Ti should be included to capture the full dynamics, 
as it is in any approximation scheme. However, the noise properties that we aim at discussing in this paper emerge 
already at the zeroth order in r^. Our scope is to design a minimal model that accounts for the non-trivial properties 
of extrinsic noise. 

As it is well known, the stochastic integral associated to the multiplicative noise term in ^ is not uniquely deflned 
when white noise is assumed [l2j . and its evaluation needs to be supplemented with a prescription on how to deflne 
the discretizcd form of The assumption that the white noise pH)) is the limit of a correlated noise when the 
correlation time goes to zero resolves the ambiguity, and leads naturally to the adoption of the so-called Stratonovich 



prescription 1^. The same prescription is also valid when an elimination of fast variables is implemented, as in the 
reduction to the Michaelis-Menten form, as long as the equilibration time of the fast variables is anyway smaller than 
the correlation time of the noise [21]. 

A useful strategy to solve (O is to derive the corresponding Fokker-Planck equation [sj], and then compute the 
stationary probability density function. The Stratonovich form of the Fokker-Planck equation for the stochastic 
process (0) with the correlation (fTU)) is 

^-T.-^A [/'(^' + P'M^^ *)} + E { [=^-5.fc(c, /i).9,.(c, /i)]p(c, t)} , (11) 
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where p{c, t) is a short notation for p{c, i|co, to)j and 



J=l k=l 



3 = 1 k=l 



M 



' d ^ 



(12) 



is known as the Stratonovich drift, here specified for the kinetics 

Under the assumption of natural boundaries the Stratonovich stationary probabihty distribution can be eval- 
uated to be: 



where 



'^''P I / E ^•'^ E (c, M) 



d 



fkic,fi) - E -Q^Dkj{c,n) 



(13) 



Efc^ifc (c, (c, /i) and /i (c, /x) = (c, /x) + G,; (c, /x) . 



(14) 



The maxima of the stationary probability distribution ps (c) are assumed to correspond to the macroscopic stable 
steady states of the pathway when noise is present, as they represent the states where the pathway is found for most 
of the time. Other indicators, such as the first moment, do not necessarily reflect well the behaviour of the system 
[2^ . For this reason we shall focus on calculating the most probable states of the pathway, under the assumption 
that they represent the proper continuation of the deterministic states when noise is switched on. By differentiating 
eq. (jl3p the extrema of the probability distribution are given by the expression: 



f(c„i,/i) - G(c,„,/2) = 0. 



(15) 



The original fluctuations on parameters do not simply appear on the concentrations as mere fluctuations around their 
zero-noise values. Rather their effect is that of flrst shifting the value of the steady state concentrations, and then 
superimposing fluctuations upon them. This mechanism originates directly from the assumption of the noise being 
external. The small correlation present (in any case smaller than the other typical time scales of the system) makes 
the noise develop correlations with the dynamical variables at the same time [7| . As a consequence the multiplicative 
noise term in eq. ([9]) develops a non-zero average. The Stratonovich prescription is a way of extracting the non-zero 
average contribution in the limit of zero, and the corresponding modiflcation of the dynamics. This modification is 
reflected in the change of the probability distribution with respect to the deterministic case, as well as of its zeroes, 
as shown by the condition In general this accounts for a change in position, number, and stability properties of 
the steady states with respect to the deterministic solution. 



3. IMPLICATIONS FOR METABOLIC CONTROL ANALYSIS 
3a. Standard Metabolic Control Analysis 

Metabolic Control Analysis (MCA) is a powerfuU framework to relate local variables, such as enzymes activities or 
rate constants, to systemic variables, such as fluxes or concentrations. It focuses on establishing how the flrst exert 
control on the latter. 

Ideally, MCA is based on a simple sensitivity analysis, which can be schematized as follows: i) Let the system relax 
to its steady state; ii) Apply a small perturbation; in) Wait for relaxation onto the new steady state; iv) Measure the 
change in global variables (fluxes, concentrations, or other systemic variables). 

The global effect of local changes is well described by a set of control coefficients, such as flux control coefficients, 



and concentration control coefficients, 
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These can be derived in the approximation of small perturbations from the kinetic description of the metabolic network, 
for any given pathway, either branched or unbranched ^ISj. Notice that the form of control cocfhcicnts presented here, 
eqs. and (fT7|) . contains the normalization factors Vj/ Ji and Vjjci. This form allows the characterizations of 

fractional changes and is suitable for the cases when none of the J^s or c^s is zero, which might happen in principle 
for certain combinations of the parameters. In that case unnormalized control coefficients may be used [isj . 

Control coefficients describe well the properties of the system as a whole. This particular fact is expressed through 
the so-called Summation Theorems, which relate coefficients for both concentrations and fluxes to each other: 

M 

Y.^l]=^^ (18) 

and 

M 

yci^^\. (19) 



Equation (fT5|) states that for any metabolite in the network the sum of the concentration control coefficients of all 
rates equals zero. Similarly, for any flux the sum of all flux control coefficients equals 1, as implied by eq. (|19p . 
The fact that these sums are constrained demonstrates that control is shared among all rates in the pathway, and 
highlights the global properties of the network. 

Local properties are equally important. In particular the response of individual rates to the perturbation of both 
internal and external variables may be very informative. This sensitivity is expressed in terms of so-called elasticity 
coefficients, whose definition is very similar to that of control coefficients: 

^v, ^ ^ ^In^'^ ^20) 

^ Vi dz 9 In z 

Here z can be a dynamical variable, such as an internal metabolic concentration, or an external parameter, such as 
a reaction constant or a fixed external concentration. When z is an internal metabolic concentration. Connectivity 
Theorems for both fluxes and concentrations can be derived [20l |: 

M 



E^^<^=0' (21) 
1=1 

M 

Y.C:in:^^-S,,. (22) 



1=1 

Here Sij is the Kronecker symbol, equal to 1 for i — j, and otherwise. Together with Summation Theorems, 
Connectivity Theorems allow for the calculation of the control coefficients from the knowledge of the elasticities, and 
therefore integrate between a local and a global description of the network. 

In the present context we are also interested in the case when z is a parameter, such as any of the fXjS introduced in 
the previous section. Perturbation of parameters is described at the systemic level in terms of response coefficients, 
such as concentration response coefficients, 

~ — 7) — ^ 7r\ — ' ^^'^) 



and flux response coefficients. 



^j, ^ Mj ^ din J, ^2^^ 



The interplay among response coefficients, control coefficients, and elasticities is formalized by the Partitioned Re- 
sponse relations, for both concentrations and fluxes. 



M M 

^S=E^:<; i?;!; ^E^^'^m'- (25) 

i=l i=l 
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Eqs. (PSj) show that the system response originates locally through the sensitivity of local rates to the external 
perturbation, and is transferee! to the systemic level through the mediation of the control that each rate has on 
the system variable of interest. Response coefficients are therefore determined by both local sensitivities to external 
effectors, and control exerted by individual rates. The set of control coefficients, response coefficients and elasticities, 
obeying Summation and Connectivity Theorems and/or Partitioned Response, represents a useful framework for 
connecting local and global properties of metabolic networks. In the next section we shall see that this mathematical 
structure is robust to extrinsic noise, but at the same time is sensitive to it through the acquired explicit dependency 
of all coefficients on the noise intensity. 



3b. Metabolic Control Analysis with extrinsic stochasticity 



One way to derive control, response, and elasticity coefficients, and their structural relationships, is to use the 
perturbative procedure mentioned in the previous section, and work out the corresponding linear approximation 
for small perturbations [l6| . However, there is another elegant and compact way of deriving such properties when 
concentrations and rates are homogeneous functions. 

Say that / is a function of n variables xi, ...a;„, which satisfies a scaling relationship, such as: 



/(A"^xi,...,A""x„) = AV(a;i, 
Then it is possible to prove that the following equality holds true: 

d\nf 



i=l 



' dlnxi 



7, 



(26) 



(27) 



and / is said to be a homogeneous function of degree 7 (Euler's Theorem). Homogeneity of concentrations and rates, 
of degree and 1 respectively, has been shown to lead to Summation and Connectivity Theorems [l^, • 

In order to follow this approach, let us consider the equation that determines the maximum probability concentra- 
tions, eq. ITT 



M 



B. N M M 



" d fdviY 




'dvh' 


dck \d^ij) 







SaSkh =0, i = 1, ■■■,N. 



1=1 j = l k=l 1 = 1 h=l 

It is not restrictive to assume specificity of the parameters with respect to the rates, namely 

dvh 



dvi , 

This allows us to get rid of the sum over h in 
kinetics of w;, with R = Ri + R2 + ■■■ + Rm , eq- 



M 



(28) 



(29) 



Also, by defining Ri as the number of parameters entering in the 
can be rewritten in the more compact form 



Siivi = 0, 



where 



and 



VI -vi, 



Ri N 

EE^ 

j=i k=i 



' d / dvi Y 




' dvi ' 


dck \d^ij) 







i = l,...,N, 



I = 1,...,M 



Sk 



(30) 



(31) 



(32) 



Equation (|30p allows for a particularly transparent interpretation. Every rate present in the deterministic system vi 
appears together with the Stratonovich drift vf , made up of as many terms as the number Ri of parameters affected by 
noise that enter the definition of vi. Notice that in the case when the dependency on concentrations and parameters 
affected by noise is linear, the drifts ((5^ have the same analytical form as the zero noise rates vi. In this case the 
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noise simply causes a renormalization of parameter values. In contrast, for general nonlinear dependencies, eq. ([5^ 
defines interaction terms not necessarily of the same form as those defining the zero noise system. 

Requiring that concentrations be homogeneous functions of degree is equivalent to require the invariance of ([501) 
under rescaling of the rates. In turn this rescaling needs to be defined consistently so that the rates themselves be 
homogeneoeus of degree 1. In this way, while individual rates undergo a scale tansformation, concentrations will stay 
invariant. 

It is clear that the scaling properties of the Stratonovich drift, eq. (|32)) . depend on the specific form of the 
deterministic rates, and on the parameter that has been perturbed by noise. We shall focus on three different cases, 
as representative of most situations in metabolic networks: i) rates obeying the law of mass action, ii) rates described 
by Michaelis-Menten kinetics, and in) rates described by reversible Michaelis-Menten kinetics. 

Consider first a rate, wiin, having the form given by the law of mass action, and let /iun be one of its parameters 
affected by noise. Since vn^ is linear in /iun, the rescaling /iun A/iun implies wun Auiin. Then the associated 
Stratonovich drift vf:^^ can be made to scale as wun only if — > Ae^ where is the intensity of the noise acting upon 
/xiin- Therefore the rate scales as A, that is it is homogeneous of degree 1, if 

v\in -> Awlin and Ae^. (33) 

Next consider rates of the irreversible Michaelis-Menten form, umm = VyaayiS/ {Km + <5'), with S the substrate 
concentration, and Vmax and Km maximal velocity and Michaelis-Menten constant respectively. The general form of 
the scaling properties of wmm is 

I'MMlA^V^nax, A°i^m) = A"wMM(Knax, K^)- (34) 

This implies wmm Awmm, when Vmax A Vmax and Km Km- To find the scale transformation of the associated 
Stratonovich drift we consider separately the two cases, when the noise acts either upon Vmax, or upon Km- If 

Vraax is affectcd by a noise having intensity ev, the corresponding w^m scales as wmm again if ey scales as A. So that 
the rate vmm scales as A if 

VMM — > A^iMM and ey — > Xsy- (35) 

Notice that ((35l) is the same set of scale transformations as (|33|) . This is not surprising since the dependency of any 
Michaelis-Menten type of rate on maximal velocities is in fact linear, and therefore the previous case must apply. 

In contrast, if Km is the parameter of vmm perturbed by a noise with intensity Ek, then the corresponding w^jm 
scales as X^Sk when Vmax AT^iax and Km — * Km- Therefore, imposing sk to scale as A^^ makes w^m scale as 
A, which is what we want in order for "Dmm to be homogeneous of degree 1: 

"MM AwMM and sk X~^£k- (36) 
The same result can be obtained also in the case of the reversible Michaelis-Menten rate, 

vrmm = -j^ — g — — , (37) 

where V^-,ax and V^^^^^ are maximal velocities for the forward and backward reactions respectively, and Ks and Kp 
are the Michaelis-Menten constants associated to the substrate S and the product P. The scaling properties of wrmm 
are similar to the irreversible case: 

WRMM(A"Vmax, A"Vmax, A"i^s, XPKp) = A"wrmM (Vjjjax: Kl?ax: ^S, Kp). (38) 

As before, if V^^ax and/or Vmax are affected by noise, urmm scales as A if 

Vrmm Aurmm and sy Aey. (39) 
On the other hand, if any of the Michaelis-Menten constants are affected by noise, then the scaling 

Vrmm Awrmm and ek A^^e^. (40) 
makes sure that urmm is homogeneous of degree 1. 
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In summary, under rescaling vi — > Aw; , the rates vi are homogeneous functions of degree 1 if the noise intensities of 
any parameter appearing hnearly in the rates are rescaled as A (including maximal velocities of possible Michaelis- 
Menten rates), and the noise intensities associated to Michaelis-Menten constants are rescaled as A~^. Also, because of 
the invariance of pop under this same rescaling, the maximum probability concentrations are homogeneous functions 
of degree 0. For other possible types of parameter dependencies, the proper scale transformation need to be sought 
case by case by applying similar arguments. 

By using Euler's Theorem, the homogeneity of concentrations and rates leads directly to the Summation Theorems 
for the respective control coefficients. In particular, from the scaling relation 

c"(ACi,...,ASAf) -c"(^i,...,«Af), (41) 

and by using ^7} . we readily obtain 

M 

d In Vj 



that is 

M 

E<-0. (43) 
i=i 

Eq. (j43p proves that the Summation Theorem for concentration control coefficients is still valid when stochasticity is 
present in the form of noise on parameters. 

Similarly, by defining the macrosopic fluxes as the fluxes evaluated in correspondence of the most probable concen- 
trations, J™ — ■D(c™), it is clear that 

J"(Ai}i,...,A{iM) = AJ"(«i,...,«m). (44) 
Euler's Theorem may be invoked again, and implies 

i=i 

which states that the Summation Theorem for flux control coefficients is valid as well, even in presence of extrinsic 
noise. 

The substitution of the zero noise rates vi with the rates pip , which obey the same scaling properties, allows also 
to derive in a standard fashion Connectivity Theorems 14 1 and Partitioned Response [l^ : 

M M 



1=1 



M M 



E 4 =0 E 4 = (46) 

1= 

and 

^S=E^r4 (47) 

i=l i=l 

It is interesting to note that new elasticity coefficients can be defined, which describes the local effect of perturbations 
of the noise intensities, 

together with new response coefficients, for both concentrations and fluxes, which in turn describe the effect of the 
noise on the sytemic variables: 

Rl^ and Bi] (49) 
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k ^ k ^ 

FIG. 1: The elementary reaction pathway chosen as an illustration of the efTect of extrinsic fluctuations. 



Summation Theorems P5)) and (|45p . Connectivity Theorems and Partitioned Response (|T7)) show that the 

mathematical structure of the network is robust against extrinsic noise. This property is non trivial, because it might 
be expected that extrinsic noise might break the homogeneity of the system. Partitioned Response originates in 
fact from the breakdown of homogeneity of fluxes, but in the same way as in the deterministic case. The fact that 
the Stratonovich drift obeys the same scale transformation as the zero noise rates guarantees that these properties 
are preserved. However, it should be noted that individual control coefficients, as well as response coefficients and 
elasticities, depend explicitly on the noise intensity. This feature demonstrates that extrinsic noise can in fact act as a 
control mechanism. It modifies the local rates, as formalized by the corresponding elasticy coefhcients, and propagates 
up to the systemic level to affect global variables. 



4. AN EXPLICIT EXAMPLE 



As an example of the theory presented, let us consider the chemical reaction depicted in Fig. [TJ Despite its 
simplicity it serves well the purpose of illustrating the notion of control by noise. 

Let us define the concentrations of the species A, B, and X as a, b, and x respectively. Let us assume that the 
constant k2 is the only parameter undergoing stochastic fluctuations. This leads to the following Stratonovich SDE 

^ = v,-v2-e'^'xat), (50) 
at 

where ^(i) is Gaussian white noise. The deterministic reactions vi and V2 are given in terms of the rate constants 
defined in Fig. (TJa as 

vi = kia — k-ix, (51) 

V2 = k2X - k^2b. (52) 

Notice that Equation ([50]) is exact at all orders in e, because of the linear dependence of the kinetics on k2- In 
general, this dependence may be nonlinear, and the equation corresponding to (|50p in the nonlinear case should be 
considered valid only at the leading order in e. 

By following Section 2, we compute the stationary concentrations and fluxes in order to derive the corresponding 
control coefficients. In one dimension, with extrinsic fluctuations acting upon only one parameter, eq. (|13p becomes: 

p,{x) = ^e^p(l r 4^dx'] , (53) 



gix) \e J g'^ix') 

with f{x) ~ vi{x) — V2{x), g{x) — ~x and TV a normalization constant that can be directly calculated by imposing 
Ps{x)dx = 1. The explicit expression of ps{x) is then 

"^("^ = [^^) r((fc_r + MA) [ ^ j ■ ^''^ 



The extrema of Ps{x), as from eq. ((T5|), result in wi — W2 — ex = 0, which gives 

fclfl + fc_2& 

k-i + k2+e 

Notice that all the moments of the distribution ([54]) can be calculated, and in particular we have: 



(55) 



kia + /c_26 

dx X ps{x) = — , e<fc_i + fc2. (56) 

k-i + K2 - e 
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FIG. 2: Left panel: Dependency over time of the concentration x, as given by eq. (|50|) . Parameters have been set arbitrarily to 
the following values: fci = 1 min~^, = 1 min~^, fc2 = 2 min~^, fc_2 = 1.2 min~^. External concentrations have been taken 
as a = 1 mM and b = 2 mM, while the noise intensity has been fixed to e = 0.5. The dashed curve corresponds to the solution 
of the deterministic equation, the spiky and dotted lines show the value of the concentration averaged respectively over 10 
and 1000 realizations of the noise. The thick straight line is the value of the average concentration has from eq. (|56|l . namely 
(x) — 1.36 mM with the chosen parameters. Right panel: The stationary probability distribution (|54|l (solid line), versus the 
estimate coming from the numerical solution of eq. (I50|l . with the parameters as above. The maximum is at Xm = 0.97 mM, 
whereas the deterministic value is Xm = 1.13 mM. 



Fig. [5] shows a direct numerical simulation of the SDE ([511)) with the Stratonovich prescription. The numerical 
integration scheme is standard (see for instance [ll[). In the left panel of Fig. [5] comparison is shown between the 
deterministic solution, namely the solution of (|50p with the noise term switched off, the full stochastic solution, and 
the expectation (x), averaged over 10 and 1000 realizations of the noise. The theoretical value of (x) as given by 
([56| is also plotted for comparison. In the right panel of the same figure we show the probability density ((54)) and its 
numerical evaluation from the solution of eq. ()50p . Notice that the maximum of the probability density has undergone 
indeed a shift with respect to the deterministic value, by changing its value from Xm = 1.13 mM (deterministic) to 
Xjn = 0.97 mM (stochastic). The same is true for the expectation (x), whose value has changed from (x) = 1.13 mM 
(deterministic) to (a:) = 1.36 mM (stochastic). 

In summary, the effect caused by extrinsic fluctuations affecting the parameter ^2 is a shift in the value of the 
maximum of the probability density, which induces a shift also on the value of its first moment. Notice that this 
mechanism corresponds effectively to resetting the value of the parameter ^2 to ^2 + e. This resetting is non-trivial 
in that it rests on the action of the fluctuations only, which shift the value of k2 effectively perceived by the system. 
Nonetheless the bare ^2 is blocked to its zero- noise value. Due to the linearity of the system, this is the only effect, 
and more subtle modifications, such as in the number of stationary states of the probability distribution, or in their 
stability properties, are not present. 

The same resetting appears in the expression for the flux as well, which becomes thereby explicitly dependent on 
the noise intensity. In fact, the rates ()HT)) read in this case: 

vi = kia — k^ix, (57) 
W2 = (^2 + e)x - fc_2&, (58) 

and their direct evaluation at the value x^ results in: 

fc_i + k2+e 

The noise dependency of the maximum probability concentration (|55)) and of the corresponding flux ()5^ can be 
interpreted within Metabolic Control Analysis. To this aim, we compute now control coefficients, response coefficients 
and elasticies, and verify that the Summation and Connectivity Theorems are fulfilled, as from Section 3b. 

The concentration control coefficients result in 



^ fci(fc2+£)Q-fc-ifc-2^ ,„„s 
'^^ (fc_i+fc2+e)(fcia + fc^26)' ^ ' 
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FIG. 3: Concentration {left panel) and flux {right panel) control coefiicients as from eqs. (|60 p -(|63 p as functions of e. Parameters 
setting is as in the caption of Fig. [S] 



and similarly for the flux J™, 



„^ _ fci(fc2 + e)a - fc_ifc_2fc 
" ~(fc_i + fc2 + £)(A:ia + fc_26)' 



J ^ k2+e 

fc_l + fc2 + £' 

fc_i + fc2 + £' 



It is immediate to verify that the Summation Theorems hold true, as expected: 







and 



Ci + C:' 



1. 



(61) 

(62) 
(63) 

(64) 



In Fig. [3] we present a numerical evaluation of expressions as functions of the noise intensity e, with 

parameters chosen as indicated in the caption of Fig. [21 As eqs. (155]) and (|59p show, the effect of the noise is to allow 
an extra tuning mechanism for modulating the systemic variables. Futhermore, noise may also allow for a change of 
the control that individual rates exert on the whole pathway. In other words, the dependency of the control coefficients 
on the noise makes it possible for the system to rebalance controls, and possibly to explore dynamical regions that 
might not be accessible given the standard parameter values. Already in this simple linear system, this effect can be 
so dramatic as to invert the sign of the control of the two rates on the concentration, as the left panel of Fig. [3] shows. 
Of course this extra tuning is always constrained by the respective Summation Theorems. 

As proven in section 3b, Connectivity Theorems are also preserved. It is a simple matter to compute the elasticities. 



fcifc_ia + fc_ifc_26 
ki{k2 + s)a — fc_iA;_2&' 
52 _ ki{k2 + e)a + fc_2(fc2 + e)b 
^ ki{k2 + £)a — k-ik-2b 



and verify that 



CI + C? K'^-l and Ci^ + C4 <^ = 0. 



On the other hand, elasticities and response coefficients with respect to noise turn out to be 



e{kia + fc_26) 



ki{k2 + s)a — fc_iA;_2^' 



(65) 

(66) 
(67) 

(68) 

(69) 
(70) 
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and 

K = -1 . (71) 

nJ ^ £fc_i(fcia + fc_2fe) .„2) 

^ (fc_i+fc2+e)(fci(fc2+£)a-fc-ifc-2&)' ^ ' 

and obey the Partitioned Response: 

i?f = C? <^ + = eg and i?'^ = Ci Trf^ + <^ = Ci np . (73) 

The elasticity (|70p highhghts the local effect of noise, and the response coefficients ([TT]) and (17^ for concentrations 
and fluxes are a manifestation of its global effect at the systemic level. While the action of noise on control coefficients 
is only a modification of the existiting ones, which become noise dependent, the elasticity (|70p and the response 
coefficients ([7T|) and ([7^ are truly "new" coefficients, which are turned on by the presence of noise. Nonetheless 
they obey Partitioned Response in the standard fashion. In this sense, the mathematical structure of the theory is 
generally preserved. 

5. CONCLUSION AND OUTLOOK 

Metabolic Control Analysis is a fundamental framework that links both conceptually and practically local and 
global properties of metabolic networks. In this paper we propose a way of generalizing MCA, and its predictive 
power at the systemic level, when fluctuations are present in control parameters, and the kinetics are affected by 
so-called extrinsic noise. 

The proposed generalization relies on the acquired explicit dependency of all standard MCA coefficients on noise 
intensity, and on the definition of new elasticity and response coefficients accounting for noise perturbations. Even 
though Summation and Connectivity Theorems, as well as Partitioned Response, are shown to be robust against noise, 
noise introduces an extra tuning mechanism, capable of modifying the local and global properties of the pathway. 

In fact, through the noise dependency of the MCA coefficients, we aim at putting forward a new way of looking at 
noise in metabolic networks. We propose to go beyond the interpretation of noise as a mere experimental nuisance, 
and to explore its potentiality to exert external functional control onto the system. As a result of our analysis, 
the behaviour of the pathway can be regulated by the insertion of ad hoc random fluctuations on properly chosen 
parameters. 

Intrinsic noise, related for instance to the low copy number of molecules or to protein conformational fluctuations, 
can also be described by this approach. From within a modelling perspective, deciding whether noise is internal or 
external to the system is so delicate an issue as defining the pathway of interest itself. In an integrated description of 
the full network all noise has an internal origin. However, the description of reduced modules, as opposed to the fully 
integrated approach, requires the adoption of effective kinetics, where intrinsic noise is translated into extrinsic noise 
on parameters. The predictions of these effective kinetics are the ones, which should be compared with experimental 
data, or to full stochastic simulations. 

This is particularly relevant when the system is characterized by a time-scale separation between slow and fast 
variables. If the latter are eliminated by standard techiques [9], their dynamical contribution is rewritten in parametric 
form. Accordingly original intrinsic fluctuations acting on fast variables are perceived as extrinsic noise in the reduced 



system. As long as the reduction procedure is carried out before the approximation of white noise 2l| , the reduced 
dynamics is modified by the Stratonovich drift. This scenario is confirmed for instance in where an improved 
version of the Linear Noise Approximation combined with the elimination of fast variables in a simple metabolic 
network shows a very good agreement with the results of a direct stochastic simulation Q . 

In all other cases intrinsic noise can be described by standard tools such as Master Equations, and possibly their 
Fokker-Planck [31] or Linear Noise approximations [91]. Effective Langevin descriptions can then be derived, but, 
apart from some exceptions [tI, these are generally characterized by zero average multiplicative noise terms, with no 
Stratonovich drift associated. This implies that no control mechanism can be expected. 

A recent work by Levine and Hwa [1^ confirms this prediction. The Authors of study the effect that intrinsic 
noise, as associated to copy number fluctuations of metabolic species, may have on metabolic pathways. The con- 
clusion of the Authors is that noise "does not propagate" , being the fluctuations of the metabolites in the pathway 
uncorrelated from each other. This is fully consistent with our framework. However possible correlations among 
metabolic fluctuations might in fact arise when fluctuations of enzymes common to two or more metabolites are 
considered. Our treatment of extrinsic noise aims to address this particular case. 
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In this respect, the possible internal origin of extrinsic noise suggests an even more interesting role of stochastic 
fluctuations. Robustness of biochemical systems to perturbations is often thought to be a fundamental property, and 
selection of robust traits has been proposed as a driving principle of evolution [l9| . Noise of course is no exception, 
and the capability of preserving functionality despite noise is certainly a fundamental feature of living systems (see 
for instance 3 for an example). A different perspective, however, emerges from considering cellular noise as a genetic 
response to chemical or physical constraints, as represented for instance by the values of reaction constants in a 
biochemical network. As we have shown, the operational point of the network is determined by the parameter values 
jointly with the intensity of the external noise acting on the system. However, as discussed above, external noise may 
have its origin in intrinsic noise, associated for instance to concentration fluctuations. It is a suggestive picture that 
biochemical networks would have evolved in such a way that concentration fluctuations allow the network to gain 
access to otherwise dynamically forbidden regions [3^ . Noise then would appear as a result of positive selection on 
gene expression levels rather than as an external factor causing "defensive" architectural evolution. To what extent 
evolution can positively select for stochastic dynamics as opposed to deterministic ones is a fundamental and intriguing 
question. 

The theory proposed in this paper applies to both linear and nonlinear dependencies of reaction rates on parameters 
and concentrations. In the linear case the analysis is particularly simple. The effect of the noise is perceived by the 
system as a renormalization of the relevant parameters. 

When nonlinear dependencies are considered, may be that the effect of the noise cannot be reabsorbed as a renor- 
malization of bare parameters only, and the overall dynamics is expected to be generally more complex. This case 
has been object of intense study, and the theory of noise-induced transition has emerged ^22]. Because of the nonlin- 
earities, the effect of the noise is not limited to a shift in the extrema of the probability distribution, to which our 
theory applies for small enough noise intensities, but new extrema may appear, and the system can jump from one 
to the other driven by random fluctuations. Even though this mechanism is well understood at the kinetic level, at 
the systemic level its description requires a proper MCA, which needs to be valid for large jumps between stationary 
states. This is the subject of some recent literature P, Hj, and the corresponding analysis when extrinsic noise is 
present will be the subject of a further study. 

As a final remark, we would like to mention also the recent extensions of MCA that have been discussed in 



27| to account for spatial dependencies. These are relevant for instance in signal transduction networks, where 



protein diffusion plays an important role in the signaling process [30[. Modelling these systems relies on defining 
proper Reaction-Diffusion dynamics, for which the standard MCA needs to be extended accordingly. Our paper 
follows conceptually the same spirit of (27| . as it aims at identifying on general grounds all processes that may have 
relevant biological implications. The simultaneous treatment of both noise and spatial dependencies in stochastic non 
homogeneous systems is in turn a combination of the theory of [27] and the theory presented in this paper. This will 
also be the subject of a further study. 
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